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Abstract. We describe a class of algebraically solvable SUSY models by 
considering the deformation of invariant polynomial flags by means of the Darboux 
transformation. The algebraic deformations corresponding to the addition of 
a bound state to a shape-invariant potential are particularly interesting. The 
polynomial flags in question are indexed by a deformation parameter m = 1, 2, . . . 
and lead to new algebraically solvable models. We illustrate these ideas by 
considering deformations of the hyperbolic Poschl- Teller potential. 



PACS numbers: 03.65.Fd, 03.65.Ge 
1. Introduction 

Our purpose in this paper is to show how new classes of exactly solvable 
supersymmetric quantum mechanical Hamiltonians arise in a natural fashion from the 
application of the Darboux transformation to classes of second order linear differential 
operators that preserve flags of vector spaces generated by univariate polynomials. 
These new Hamiltonians and their bound states have closed analytic expressions in 
terms of elementary functions, and their qualitative behavior is both natural and 
significant from a physical point of view. 

Recall that the Darboux transformation starts from the knowledge of a formal 
eigenfunction of a Schrodinger operator, which is used to factorize it as a product of 
first-order operators. Depending on whether this formal eigenfunction or its reciprocal 
are square integrable, one obtains forward or backward Darboux transformations, in 
which the supersymmetric partner Hamiltonian is obtained by reversing the order of 
these factors. The principle of our approach is to consider only those factorizations for 
which the effect of the Darboux transformation on functions is to map polynomials to 
polynomials. These are given a simple characterization in our paper in terms of the 
starting formal eigenfunction. We will refer to this transformation as the algebraic 
Darboux transformation. The new Hamiltonians obtained in this fashion will be 
exactly solvable in the precise algebraic sense that they will also admit complete 
invariant flags of polynomial subspaces. 

When a parametrized family of potentials is closed with respect to the forward 
Darboux transformation, it is said to be shape-invariant [1] and in this case the 
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iteration of the transformation furnishes a complete description of the spectrum and 
eigcnfunctions. For the shape invariant potentials, the underlying invariant flag of 
polynomials is the full polynomial module. Thus, to obtain deformations of a shape- 
invariant potential one must consider the two-parameter family of backwards Darboux 
transformations. These were first applied to the harmonic oscillator in [2], while the 
general theory was developed in [3,4]. However, as noted in [6], the general form of the 
deformed potential can only be expressed by a formal power series, or as the integral 
of eigcnfunctions of the original Hamiltonian - in contrast to the original potential, 
which is an elementary function, with bound states also described by elementary 
functions. Our main emphasis in this paper is to obtain examples of exactly solvable 
potentials which lie outside the shape-invariant class, which can be expressed in closed 
analytic form, and which have qualitative properties that make them relevant to the 
description of physically realistic situations. These are obtained by the application of 
the backward Darboux transformation to shape invariant potentials corresponding to 
special values of the parameters. There is a countable infinity of algebraic backwards 
transformations of a shape- invariant potential [7], indexed by an integer m. We will 
show the precise manner in which the m'^ algebraic backward transformation deforms 
the invariant polynomial fiag of a shape invariant potential, and we calculate the 
explicit basis of the deformed flags for the cases m = 1 and m = 2. As an illustrative 
example we discuss in detail the algebraic deformation of the hyperbolic Poschl- Teller 
potential. 



2. Darboux trEinsformations 

2.1. The self-adjoint case. 
Consider the Schrodinger operator 

H = -dxx + u, (1) 

where u{x), x €.Ris continuous, real- valued and bounded from below. Consequently, 
the restriction of _ff to a certain dense subspace D{H) C L^(]R) is a self-adjoint 
operator. Consider a formal eigenfunction 

H[4>] = Ao(/). 

The key idea of the supersymmetric or Darboux transformation is the fact that to 
every ^ there corresponds a factorization of H as 

H-\q = a* a, (2) 

where 

A = d^- (log 4>U A* = -d, - (log 4.),. (3) 

We shall refer to 4> as the factorization function, and to Ao as the factorization energy. 
The supersymmetric partner potential is the operator defined by the commutation of 
the factors 

H = -da:x +U = AA* +Xo, U = U- 2(l0g (l))xx- (4) 

The transformed potential u is continuous if and only if </) is non-vanishing, which we 
assume from here on. In this way, H is self-adjoint and semi-bounded on some dense 
domain D{H). The operators H and H satisfy the intertwining relation 

AH = HA, (5) 
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which imphes the following relation between the eigenfunctions of the two operators: 

H[^] = X^P, H['ip] =x4j, i} = A[^]. (6) 

The spectral properties of this transformation are governed by one of the following 
three possibilities [4,5]. 

(i) Forward transformation: 4> is square integrable (and since it is nodeless, it 
must be the ground-state wavefunction of H). The operator A maps D{H) onto 
D{H), with a 1-dimensional kernel. The n"^ bound state of H is mapped to the 
{n — bound state of H. Correspondingly, the transformed spectrum differs 
from the spectrum of H by the removal of Aq, the lowest eigenvalue. 

(ii) Backvifard transformation: is square integrable. The operator A maps the 
n^^ bound state of H to the {n+ bound state of H. It is one-to-one on D{H), 
but not onto D{H)] the new ground state is not in the image. The spectrum of H 
differs from that of H by the addition of a lowest eigenvalue, namely A, with the 
ground state given by A 2-parameter family (energy and shape parameter) 
of backward transformations exist for every A strictly smaller than the infimum 
of the spectrum of H. 

(iii) Isospectral transformation: neither (p nor i/)"^ are square integrable. The 
operator A defines a linear isomorphism from D{H) to D{H). It transforms 
the n'^ bound state of H to the n"^ bound state of H. Two isospectral Darboux 
transformations exist for every A strictly smaller than the infimum of the spectrum 
oiH. 



2. 2. The general form and covariance. 

In this paragraph we will consider Darboux transformations of an arbitrary second 
order operator, the general form of which is 

T=pd,,+qd,+r, (7) 

where we assume that p{z) < on the domain of interest. The above operator is 
related to a Schrodinger operator (QJ by the change of variables 

{-p)-^dz, (8) 

and gauge transformation 

H = ePTe-P, p ^ j \p-\q - \p,) dz. (9) 

The relation to the potential is given by 

u = \pzz - \qz - \p-^ {q - ^Pz)iq - ^Pz) + r (10) 

Since gauge transformations and changes of variable are homoniorphisms of the ring of 
differential operators, the Darboux transformation is covariant with respect to these 
operations. 

Although it is customary to work in the Schrodinger gauge as in Section |2. II the 
Darboux transformation can be defined relative to a general coordinate and choice of 
gauge, as shown below. Indeed, let 

T[0] = Ao0, (11) 
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be a factorization eigenfunction. Writing 

(log0). = (12) 
we have the following factorization of T: 

T^BA + Xo, (13) 

where 

A = bd,-a = b{d,- {log cj,),), (14) 
i? = 1(^5. + + ^) - {pd, +P(log0). + q)b-\ (15) 

We define the partner operator to be 

f = AB + Ao, (16) 
and observe that the following intertwining relation holds 

TA AT. (17) 

A particular case occurs when q — 1/2 pz, i.e. when the operator is in the self-adjoint 
gauge. Taking 

6=(-p)i/2, a=(-p)i/2(log0),, 

we have B = A*, and consequently (|13() is equivalent, after a change of variables, to the 
self-adjoint factorization We also note that the mapping T i-^ T is not canonical, 
but rather covariant with respect to gauge transformations, since two different choices 
of the denominator in (|12|l 

(log0), =a/6-a75', 
will lead to partner operators T and T' which are related by a gauge transformation: 
b-'fb-Xo = {b')-' T' b'-Xo = id, - (log <\>)z){pdz + q + p(log (/.),). (18) 
To effect an inverse transformation, we factorize T with 

(f)^bexp(^~J (j^ + ljd^^ (19) 
as the factorization function. A simple calculation shows that 

f [0] = Xo4>, 

and T has the following form 

f^pdzz + qdz+f, 



where 



q = q + Pz-2p{\ogb)z, (20) 

r = - p(log b)zz + p(log b)z^ - {pz + q) (log b)z (21) 

2pa^ , , a 

- -yT + (f^- - 2<?)- +qz + 2Ao - r. 



Thus, taking 



a = p(&z — a) — qb, b = pb, 
we have (log 0)2 — a/b, and from (fn|) (fT^ we obtain 
Bi + Ao, 
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where 




(22) 
(23) 

It follows that 

f = AB + \a = b^Tb-^. 
2.3. Algebraic factorizations. 

We will say that a second-order differential operator is exactly solvable by polynomials 
(P.E.S.) if it is equivalent, by a change of variable and a gauge transformation, to a 
second-order operator T that preserves an infinite flag of finitc-dimcnsional polynomial 
subspaces 

Ml C M2 C M3 C . . . C M = UnMn- (24) 
As part of this definition we include the following assumptions: 

(El) There is a fixed polynomial codimension, which we will call to. To be more 
precise, we assume that each Ain is an n-dimensional subspace of 

Vn+m-l^{l,Z,z\...,z"+'^-'). 

(E2) There is no spectral degeneracy. The action of T is upper-triangular relative 
to a basis adapted to the above flag, and hence possesses an infinite list of 
eigenpolynomials. We assume that the corresponding eigenvalues are distinct. 

This definition is similar to the definition of exact solvability introduced in [8]. Of 
particular interest is the subclass of P.E.S. operators for which Adi = R. If this 
condition holds, we will say that the operator satisfies the algebraic ground state 
condition. 

We will now consider the following question: suppose that T is a P.E.S. operator 
with invariant flag (|24|l . what are the conditions on a factorization function (jj such 
that the partner operator T is also P.E.S. ? To answer this question, we define a 
factorization eigenfunction Hll|) to be of algebraic type whenever 



is a rational function, where without loss of generality the polynomials a = a{z) and 
b = b{z) are assumed to be relatively prime. This definition is motivated by the 
following observation. 

Proposition 1 The factorization function (j) is of algebraic type, if and only if the 
operator 

A = bdz - a = b{dz - {logcj))^), 
transforms polynomials into polynomials. 
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We would like to focus on those Darboux transformations that preserve the P.E.S. 
character, from here on called algebraic Darboux transformations. According to 
Proposition ^ these are precisely the ones in which the factorization function (j) is of 
algebraic type. Let us describe the invariant polynomial flag of the partner operator 
in each of the three cases discussed in Section ITTl 

(Al) Algebraic forward transformation: In this case AAi = {(j)), where (p is the 
factorization function. The operator 

A = (j>dz- <j)z 

deletes the ground state, and therefore the invariant flag of the partner is 

Mn = A[Mn+l]- 

(A2) Algebraic backward transformation: In this case, the new ground state (j), 
given by (|19|l . is a rational function. Writing ip = d/b, we see that the first 
condition is true if and only if there exist polynomials d = d{z), and b = b(z) 
such that 

bz q a dz bz ( d\ 
b p b a b \ bj^ 
The partner flag is given by 

Mi^ {d), Mn+i = bA[Mn] © (a). 
(A3) Isospectral transformation: In this case both A and B are linear 
isomorphisms, and we have 

Mn=A[Mn]- 

Of particular interest is the case where the partner operators represent an 
algebraic forward/backward transformation pair, and where both operators satisfy 
the algebraic ground state condition. 

Proposition 2 Let T, T, (j), a, b be as above, and suppose that T satisfies the algebraic 
ground state condition. The following are equivalent: 

(i) T ^ T is an algebraic backward transformation, with T satisfying the algebraic 
ground state condition. 

(ii) poz + {r- Ao)6 = 0. 

(iii) (j) is a constant. 

Proof. Without loss of generality, T annihilates the constants, and it is therefore of 
the form 

T = pdzz + qdz. (26) 
The implication (iii) ^ (i) follows directly from the definitions. In order to prove the 
converse, suppose that (i) holds. By ((T3|l ifTOjl we have then 

Hence, 




pdzz + qdz-^{^] (27) 
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Note that A[\\ — ~a, and therefore the first two invariant subspaces of the partner 
operator AB are given by 

{(f>) ^ Ml and {(t>,a)=M2- 

Note that T satisfies the algebraic ground state condition if and only if a/(j) is a 
polynomial. Also, by (|25|l H27|l we have 

i) 

Since (A2) implies that Aq 0, divides both a and 6, and therefore (f) must be a 
constant. 

Finally, let us show that (iii) is equivalent to (ii). By (|19|l . (iii) is true if and only 

if 

bz — a q , ^ 
= 0. 28) 

By assumption, 

= p4>zz + q^z = Ao0, and 0^ = ^(/), 



and therefore 



„2 abz\ /aN2 



or equivalently, 

- Aofe = pa 



h p ^ 

so that (ii) holds if and only if H28|) does. □ 



3. Algebraic deformations of shape invariant potentials. 

3.1. Shape invariance. 

Let us recall that a parameterized potential is called shape-invariant [1] if the forward 
Darboux transformation preserves the form of the potential while altering the value 
of the parameters. In the preceding section, we pointed out that the Darboux 
transformation is covariant with respect to arbitrary changes of gauge and variable. 
As a consequence, the notion of shape-invariance makes perfect sense for general 
second-order operators, and not just for operators in Schrodinger form. Thus, we will 
adapt the usual definition and say that a parameterized family of P.E.S. operators 
is shape-invariant if that family is closed with respect to the forward (Al) Darboux 
transformation. 

We now describe an important class of shape- invariant, P.E.S. operators. We 
define the standard polynomial flag to be: 

R = VoCPiClV2Cl ...cVnC..., P„ = (29) 

The general form of a second-order operator T that preserves the standard flag is 

T = pdzz + qdz + r, (30) 

where p = p{z) and q = q{z) are, respectively, second and first degree polynomials, and 
where r is a constant. The family of operators described by H30|l is shape invariant in 
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the above sense. The ground state is given by (/) = 1 with Aq = r, and the factorization 
is simply 

T ={pd,+q)d,+\o. 
The partner operator 

f = d:,{pdz + (?) + Ao = pdzz + qdz + r, 
retains the form H30|l . with 

q^Pz+q, f = qz+r. 

The corresponding non-singular potential forms — see © © (|10|l for the 
transformation formulas and [7] for their derivation — are shown in Tabled These 
are the well-known, shape-invariant potential families: the harmonic oscillator (I), 
the Morse potential (II), and the hyperbolic Poschl- Teller potentials (III). Since 
potentials (I) and (III) are even functions, the corresponding eigenfunctions have a 
well-defined parity. Consequently, these potentials possess two algebraic sectors, i.e. 
they are exactly solvable by polynomials in two distinct ways (see [10] for an algebraic 
explanation of potentials with multiple algebraic sectors) . The even sector corresponds 
to an even gauge factor, and the odd sector to an odd gauge factor. The parity of the 
algebraic sectors is reversed by a Darboux transformation. 





Ic 
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II 
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9 
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9 

X 


X2 


- (yl+ i)e-^ 


\{\-A^)s^c\?{l) 


i(i-A2)sech2(f) 






Table 1. 


Shape-invariant potentials 


on the Une 





3. 2. Deformations of the standard flag. 

Let a = a{z),b = b{z) be relatively prime polymomials, and let g = g{z) be a 
polynomial that divides a^, b and bz — a. Consider the differential operators 

B = g~^dz, A = bdz-a, (31) 

and note that, by assumption, the second order operator 

T = BA = pdzz + qdz - az ff"' (32) 

where 

p = bg~'^, q={bz~a)g~^, 

has polynomial coefficients. We will say that A, B constitute a deformation pair of 
order m — deg(f;) if T leaves invariant M.n — Vn-i for all ri, i.e. if T leaves invariant 
the standard polynomial flag. Deformation pairs are of interest because they provide 
non-trivial examples of exactly solvable Hamiltonians outside the shape- invariant class. 
As usual, we define a partner operator 

f ^ AB ^ pdzz + qdz. q^ -{bgz + a)g-\ (33) 
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and a partner flag 

M = TWi c 7W2 C 7W3 C . . . , Mn = A[Vn-2] © K, 

which we wifl refer to as a deformation of the standard polynomial flag. 

Proposition 3 Every M.^ is a codimension m subspace of Vn+m-i- 

Proof. By assumption, azQ^^ is constant, and dcg[bg^^) < 2. This implies that 
deg(a) = m + 1 and that deg{b) < m + 2. Therefore, Mn C Vn+m-i- Since a,b are 
relatively prime, A does not annihilate any polynomial, and hence dimA4„ = n. □ 

Proposition 4 The partner operator T is P.E.S. 

Proof. The operator T preserves the deformed flag because of the intertwining relation 
TA = AT, 

and because T annihilates M = A^i by construction. Condition (El) is true by the 
preceding proposition. We noted above that A and hence that T does not annihilate 
any polynomial. Hence is not an eigenvalue of T , which proves condition (E2). □ 
Let us also note that a deformation pair satisfies the conditions of Proposition[51 
and in particular T ^ T is corresponds to an algebraic backward Darboux 
transformation. 

Proposition 5 The deformed subspaces Mn can be characterized as the subspace of 
'Pn+m-i consisting of all polynomials f = f{z) such that g divides fz- 

Proof. First, note that for 

f = A[h], heMn-l^Vn-2, f eMn 

we have that g~^fz — T[h], which proves that fz is divisible by g. In order to prove 
the converse, let us note that the subspace of all / G Pn+m~i such that fz is divisible 
by g is rt-dimensional. However, dimA^„ = rt by Proposition which proves the 
claim. □ 
We will now show an explicit basis of = n in the cases m = 1 and m = 2. In 
the first instance, since the subspaces Vn are invariant with respect to translations, 
we may without loss of generality assume that g{z) = z. By H32|) necessarily, b = pz, 
where p = P2Z^ + Piz + po is a polynomial of degree 2 or less. In order for g to divide 
Uz and bz — a we must have 

A = {p2Z^ + piz + pa)zdz - {a2Z^ + Po), (34) 

where a2,pojPijP2 are arbitrary real numbers. Let us also assume that 02 ^ and that 
0-2 /P2 is not a positive integer. If these generic conditions hold, then the subspaces of 
the partner flag are given by 

A[P„„2] © K = Mn = (1, z\..., Z") 

The above monomial-generated subspace is exceptional in that it admits a seven 
dimensional vector space of second-order operators that preserve it [11], and 
consequently can be used to construct novel instances of exactly solvable and quasi- 
exactly solvable operators [7,10]. 

Turning to the case m = 2, we limit our discussion to the generic case of g{z) with 
distinct roots. By scaling and translating z, as necessary, we may assume, without 
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loss of generality, that g = — 1. By H32|l . in order for g to divide and bz — a we 
must have 

A = {p2z'^+piz + po){z^ - 1)5, + {p2+Pq){z^ - 3z) - 2pi, (35) 

where po,pi,p2 are arbitrary real numbers. Let us also assume that P2 + Pa ^ 
and that —po/p2 is not a positive integer. If these generic conditions hold, then the 
subspaces of the partner flag are given by 

A[7'„_2] © M = 7W„ = (I,7r3(z),7r4(z),...7r„+i(z)), 

where 

n2k+i{z) - z^'^+i - {2k + l)z, n2k{z) = z^'' - kz\ (36) 

The above polynomials tt = 7r(z) have the property that tt, is divisible by z'^ — 1. 
The resulting polynomial subspaces M.n are preserved by the following second-order 
operators: 

T3 = z'^dzz + [{1 - n)z' - 5 + n - 5., 
- (z2 - 1)9,, - 2za„ 

Ti = za,, - 2 (^1 + a„ 

5.5. Algebraic deformations of the hyperbolic Poschl- Teller potential 

The hyperbolic Poschl- Teller potential [14], which includes the class of reflectionless 
1-soliton potentials [15], has the form 

t/pT(x) = i(i-a') sech^d). (37) 

The general solution [16, Sec. 2.9] of the corresponding Schrodinger equation 

Hpt{4') = -(j)xx + c/pT = -k^4' 

can be given as 

(^PT(a:;fc,Co,Ci) = cosh(f )^-"{Co 2Fi(-f + i + fc, -f + i - fc, i; - sinh^d )) 
+ Ci sinh(f ) 2Fi(-f I + fc, -I -f I _ fc^ 3. _ sinh^d))}, 

where 2^1(1,^, c;z) also denotes the analytic continuation of the hypergeometric 
function to Re(z) < 0. For a > 1/2 , the potential (|37(l has \a — \~\ bound states 

V'pt,i(a;), Q <i < a ~ \. 

The even bound states are given by [7] 

^Pt,2j(a:)cx0p,(x;f -j-i,l,0) (38) 



oc cosh(f)^-"pj ^' "^(coshx). 



The odd ones are given by 



i^,,,23 + l{x) cc (j)^^{x- f - J - |, 0, 1) (39) 
oc sinh(|) cosh(|)5"" pj^'""^ (coshx). 
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where Pj"'''''(z) are the Jacobi polynomials. We focus on deformations of potentials 
with bound states only, i.e. we must take a > ^- I^i order to have a well defined 
backwards Darboux transform of the hyperbolic Poschl- Teller potential, we must 
consider the solutions 0pt which correspond to an energy below the spectral minimum 
and are nowhere vanishing. Since the spectral minimum is — — a)^ we take 
\k\ > ^ — a. Now, the two-parameter family of backward Darboux transformations is 
given by the transformation functions: 

6 (.-kit) |i|<2£ii±A^lM±i±ll (40) 

with the extreme values of the shape parameter t corresponding to an isospectral 
transformation [7]. In general, for an arbitrary value of the energy parameter k the 
partner potential will be defined formally by a power series. However, for specific 
values of k the log-derivative of (/)pT will be a polynomial in cosh x thus giving rise to 
an algebraic deformation. It can be shown that these algebraic deformations occur 
precisely for the following countable subset of l|4(J|) 

(l)i"T\x) = 0PT(a;, - i - m, 1,0) 

(X cosh(|)5+"pi"^'"\cosh(x)). 
The resulting deformed potentials, as given by I^J, have the form 

U^t\x) = -i(a+i)(a-H|)sech2(f)-2 (log Pi^^'^^coshx)) ,(41) 

and have been plotted in Figure 




Figure 1. Algebraic deformations Up^ (x) of the hyperbolic Poschl- Teller 
potential I41i . with a = 4 and m, = 0, 1, 2 and 3. 



More specifically, the first (m = 1) and second (to — 2) deformations have the 
following forms: 

[/«(.) = f/(°)(.) + ^-l(^, 

21 Z{ 
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where 

zi = \{{2a + 3) cosha; - (2a + 1)), 

and 



2 



2a 



(^1 - 1)^ 

where 

Z2 = i/3((2a + 7)cosha;- (2a + l)), /? = „ o 

* V oa + 6 

The algebraic backwards Darboux transformation corresponds to the first order 
operator 

p(5."+l), 

c^a; — ^(m + a+ i) sinh a; 



/^^-r'^(coshx) 



pi '-"^(coshx) 

Both the undeformed and the deformed potentials are even functions, and 
consequently the corresponding Hamiltonians leave invariant the spaces of odd and 
even functions. The Darboux transformation changes parity. In particular, the 
deformed even sector is the A-image of the undeformed odd sector. 

We now determine explicitly a basis of the invariant flag corresponding to the 
even sector of the first and second deformation. To do so, we switch to the algebraic 
variable and perform a change of gauge so that the undeformed, odd algebraic sector 
is isomorphic to the standard polynomial flag (c.f. Case IIIq of Table QJ. 

T = z{l- z)d,, + ((a - |)z + 1 - a)d, - (f - |)' , (42) 

T = e-PH^^eP, e" = sinh(f ) cosh(f )^-", 

z — cosh^(|) — i(cosha; + 1). 

In the algebraic gauge, the factorization functions for the backward transformations 
are given by 

= = z"(z - l)-3Fi"^'"^(2z - 1) 

This factorization function is of algebraic type with 

n .^ " 5 n ,Piir'^(2z-l) a 

log(A , = - + + i + a + m "] = T, 

^ 1"^ Pi"^'"^(2z-1) ^ 

a = ((i - a)z + a)p}n^'°'\2z - 1) + (1 + a + m)z(l - z)P,^jl"^^\2z - 1), 

b = z{l- z)p!n^'"\2z-l). 

A direct calculation shows that the above functions satisfy condition (|28|l . and 
therefore 0=1. The operators 

A = bd,-a, B^g-^d,, 

where 

g = Pir^'"' (2z-l), 
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constitute a deformation pair of order m. Let us consider the cases m = 1 and rri — 2 
in more detail. For m — 1 we have 

_ zi(l - zi)(zi + 2(a + 1)) „ , (2a + l)^^ _ 4(a + 1) 



2 a 

— Oz 



2(2a + 3) ' 4(2a + 3) 

Zl 4(a+ 1) ^ 2a + 1^ 



T=z(l-zR,+ (^--^2c, + 3)zi ' 2a + 3 

Zl = 2Pi^"^'"'(2z - 1) = (2a + 3)z - 2a - 2 

The operator A = b{zi)dzi — a{zi), relative to the zi variable, is of the form (|34|1 . and 
hence, the partner operator T is P.E.S. with invariant subspaces 

yw„+i = {i,z^,zl,...,z^). 

For m — 2 we set 

Z2 = i/3 ((2a + 7)z- (2a + 4)), 

so that 



(2z-l) 



3(« + 2) 
2(2a + 7) ^ ^ 



1) 



Consequently, 



A: 



B 



-6/?z| 3(2a+l)z2 



2a + 5 



3a + 6 



3f3 



3(a + 2)2(z|-l) 



2(2a + 7)2 



+ 



3(a + 2) / (2a + 3) 3 



(2a + 7)2 
2(2a + 7) 



2/3 



{zi - 3z2) - 2a - 1 , 



Ha + 2){zi 



T = z{l- Z)dzz + 



1)^- 
3/3(a + 2) 
2a + 5 



Z2 



4(2a + 3) Z2 
2a + 7 - 1 



2(2a + l)(z| + i) 



(2a + 7)(z22 



dz 



In the same manner, it can be seen from (|35|) that the partner operator T is P.E.S. 
with invariant subspaces (c.f. 



Mr 



{l,TT3,{z2),TTi(z2), ■ . . 7r„+l (^2)) • 



4. Discussion 



In this paper we have analyzed the connection between the Darboux transformations 
and exact solvability by polynomials, i.e. the fact that a certain Hamiltonian operator 
after a change of variables and a gauge transformation admits an infinite flag of 
invariant polynomial subspaces. Since operator composition and factorization are 
covariant with respect to changes of variables and gauge transformations, concepts like 
the Darboux transformation or shape-invariance are also covariant. It is customary 
for physical applications to work in the Schrodinger gauge and the physicial variable 
X, but for the purposes of analyzing invariant flags it is more convenient to work in 
the algebraic variable z and in the algebraic gauge, the one in which the operator has 
polynomial eigenfunctions. 
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A general backwards transformation on a shape invariant potential will lead to 
a transformed potential whose eigenfunctions are no longer elementary functions. 
We have discussed the special class of algebraic Darboux transformations of shape- 
invariant potentials, i.e. those that preserve the exact solvability by polynomials, 
showing also how the polynomial flag is deformed by the action of the Darboux 
transformation. 

In this paper we placed our emphasis on the action of the Darboux transformation 
on the invariant flag of subspaces, rather than on the potential. In fact, many different 
potentials have the same invariant flag, e.g. all the shape-invariant potentials preserve 
the standard polynomial flag. We have analyzed the deformations of the Poschl- Teller 
potential, but similar deformations exist for other shape-invariant forms. 
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